Statistical mechanics of a Feshbach coupled Bose-Fermi gas in an optical lattice 
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We consider an atomic Fermi gas confined in a uniform optical lattice potential, where the atoms can pair into 
molecules via a magnetic field controlled narrow Feshbach resonance. The phase diagram of the resulting atom- 
molecule mixture in chemical and thermal equilibrium is determined numerically in the absence of interactions 
under the constraint of particle conservation. In the limiting cases of vanishing or large lattice depth we derive 
simple analytical results for important thermodynamic quantities. One such quantity is the dissociation energy, 
defined as the detuning of the molecular energy spectrum with respect to the atomic one for which half of the 
atoms have been converted into dimers. Importantly we find that the dissociation energy has a non-monotonic 
dependence on lattice depth. 



I. INTRODUCTION 

The use of Feshbach resonances with uhracold atoms in op- 
tical lattices provides a new avenue for creating molecules. 
This system has seen experimental realizations in a variety of 
systems with bosonic and fermionic species U |^ [H 0, Et]. 
For the case of fermionic atoms the Feshbach assisted conver- 
sion is to molecular bosons with strikingly different behavior 
in the degenerate regime. In this article we study the statis- 
tical mechanics of such a degenerate gas of fermionic atoms 
confined in an optical lattice and subject to a Feshbach reso- 
nance. Our interest is the chemical equilibrium between the 
fermionic atoms and their dimerized, bosonic counterparts. 

For an ultracold atomic gas a Feshbach resonance is ob- 
served when a closed channel bound state is coupled to the 
scattering continuum of an energetically open channel. The 
position of this resonance is tunable, since the different chan- 
nels correspond to different combinations of internal atomic 
states and hence experience different Zeeman shifts in an ap- 
plied magnetic field. As the closed channel bound state is 
tuned from above to below the open channel threshold, the 
resonance becomes a true molecular bound state of the two- 
body system The parameter controlling the position of the 
resonance as the magnetic field is varied is the resonance en- 
ergy i?rcsi which is the detuning of the closed channel bound 
state from the open channel threshold. 

Previously the thermodynamics of an atomic Fermi gas 
with Feshbach resonant atom-molecule conversion has been 
studied both in free space |0] and in a harmonic confining po- 
tential i^l^. For E-rcs < the dimers are stable against disso- 
ciative decay and hence are real molecules. For E^-cs > the 
dimers have a finite lifetime, but we can still discuss chem- 
ical equilibration of the gas in the steady state limit, where 
a detailed balance is established between molecule formation 
and decay. In this respect we also use the term molecules for 
the unstable dimer states, which in an ensemble will have a 
finite occupation at any instant. Note that this implies a nar- 
row Feshbach resonance, for which the scattering resonance 
is a consequence of a l ong -lived quasi-stationary state embed- 
ded in the continuum lUlll . In the structured continuum of the 



optical lattice the molecular states may also be stable against 
dissociative decay at energies above the continuum threshold, 
if their energies lie in the band gaps lfl2i[l3i[l4l[l5ll . 

For a narrow Feshbach resonance it has been established 
that the existence of a small parameter facilitates a quanta- 
tively correct description of the many-body physics based on 
a pertubative expansion iflolfTTll . Building on this insight, we 
consider for simplicity an ideal gas mixture, such that the only 
effect of the Feshbach resonance is to maintain the chemical 
equilibrium between the two species. Even this simple model 
captures the essential physics of atom-molecule conversions 
in experiments lfl6ll . and it is exact in the limit where the res- 
onance is infinitey narrow ifTTll . With a finite atom-molecule 
coupling an effective atom-atom interaction arises upon the 
elimination of the molecular degrees of freedom. We empha- 
size that even though this effective interaction may diverge 
as the magnetic field is varied across the resonance, all ther- 
modynamic quantities such as the chemical potential and the 
condensate fraction remain well behaved (see e.g. iflTi O). 
In particular, for a narrow Feshbach resonance the ideal gas 
thermodynamics remains qualitatively correct when the atom- 
molecule coupling is small but finite [fill]. 

We believe that the optical lattice introduces considerable 
new physics, and that the best way to understand the result- 
ing changes in the thermodynamics is to consider the ideal 
case first. Thus, for a narrow Feshbach resonance we an- 
ticipate that the inclusion of interactions will only impact 
our results quantitatively. Notwithstanding, for a broad Fes- 
hbach resonance interactions play a crucial role, and hence 
in that case the thermodynamics of the BCS-BEC crossover 
requires more advanced modelling than presented here, i.e. 
a full many-body theory for resonantly interacting Fermi 
atoms 111, I23, El mill IH- For a deep lattice this has 
been studied at zero temperature by Koetsier et al. i25ll . How- 
ever, even for a broad resonance a thermodynamic description 
based only on the molecular binding energy, the temperature 
and the total number of atoms can give a molecule fraction 
which agrees well with experimental data across the Feshbach 
resonance ll26ll . 

The thermal properties of ultracold atoms in optical lattices 
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have received considerable recent attention. In the degenerate 
regime ideal Bose and Fermi gases have been investigated in 



the uniform lattice 11271 12811 and in the presence of addition- 
al external harmonic confinement ll29l l30l [3lll . Interactions 
play an important role in deep lattices, and more recent work 
has examined the effect these interactions have on the thermal 
excitation generated during the preparation of bosonic Mott- 
insulating states |32, 33], and on the feasibility of achieving 
the fermionic Neel state i34ll . Finite temperature mixtures of 
atomic Bose and Fermi gases in lattices have been studied llssll 
in an attempt to explain recent experiments 1361] . 

Our system has several parameters which can be varied in- 
dependently in a numerical calculation. We only show results 
for half filling which is of the greatest interest in relation to 
current experiments. However, we have found that systems 
with filling fractions less than or equal to unity have much the 
same behavior For filling fraction larger than unity the ex- 
cited bands play a larger role, and some of the conclusions 
presented here have to be amended. For simplicity we re- 
strict our analysis to the case of a spin-balanced Fermi gas. 
There are many interesting effects associated with spin polar- 
ization, but these do not play a major role in the transition 
between atoms and molecules which we intend to study here. 
Instead the essential quantities governing the phase diagram 
are the resonance energy, i.e. the energy offset between atoms 
and molecules, the lattice depth, and the temperature. Tun- 
ing the resonance energy at a fixed lattice depth shifts the en- 
ergy spectra of the atoms and the molecules with respect to 
each other, thereby moving the point of chemical equilibrium. 
Since the atoms and molecules experience different lattice po- 
tentials, the depth of the periodic potential determines where 
the transition between the two species occurs. 



n. FORMALISM 

We consider a simple cubic optical lattice with sites 
containing A^tot identical fermionic atoms of mass with 
two internal states labeled |t) and | j), which we will refer 
to as spin-up and spin-down respectively. In this paper we 
restrict our attention to the case of an equal number of atoms 
in the two internal states 



(1) 



and we define the filling fraction ?] to be the average number 
of each kind of atom per site 



A/3 M3 2A/3 ■ 



(2) 



With the Feshbach resonance a spin-up and a spin-down 
atom can couple to a bosonic dimer with mass m„i = 2ma- 
In the following we determine the phase diagram of the sys- 
tem under the assumption that the only interaction between the 
atoms is the Feshbach resonance, which maintains chemical 
equilibrium between unbound atoms and diatomic molecules. 
Since we are assuming thermodynamic equilibrium, it is then 
meaningful to define the number of molecules iV,„ and the 



number of atoms Na as the average number of molecules and 
atoms, respectively, in a long time period where all the exter- 
nal parameters are kept constant. At all times particle conser- 
vation leads to 



Na + 2Nrn = + + 2N„ 



(3) 



where iVj and are the numbers of unbound spin-up and 
spin-down atoms respectively. 

In thermal equilibrium the temperatures of the atoms and 
the molecules are the same, T, while the condition of chemical 
equilibrium can be expressed from the atomic and molecular 
chemical potentials as 



(4) 



Since, in the spin-balanced case, the spin-up and spin-down 
atoms have identical thermal behavior, we have /^J = = /i 
and therefore /i„i = 2^. 

Our choice of energy convention is to measure all single 
particle states from the lattice potential zero of energy, and 
furthermore, we subtract the magnetic field shift E^cs from 
the molecular states. Thus the distribution functions take the 
form 

fa{E) = ,p_.,wLt^ , (5a) 



ME) 



^{E-f,)/ksT I' 
1 

^{E+E,,,~2fj,)/kBT _ I ■ 



(5b) 



As we discuss in further detail below E^ and E!^ are the avail- 
able single particle energy states for atoms and molecules, re- 
spectively, where r = 0, 1, 2, . . ., with r ~ corresponding 
to the ground state. Importantly, the molecular chemical po- 
tential is bounded from above by the lowest possible single- 
molecule energy 



Mr, 



<E 



Ere 



(6) 



since the molecules are governed by Bose-Einstein statistics. 
The density of states for atoms and molecules are given by 



Pa{E) = 2Y,S{E^E:), 

oo 

P„.{E)=Y,S{E-E:^), 



(7a) 



(7b) 



where the factor of two in (jTaJi arises because there are two 
types of atoms with identical energy spectra. The total number 
of atoms and molecules for given values of E^cs, M and T can 
be found by integrating the density of states weighted by the 
occupation over all energies. Because the energy levels are 
discrete, this becomes a sum 

Na^ f dEpa{E)fa{E) = 2^, faiEl) , (8a) 

N^^ fdE Prn{E) frn{E) = ^ • (8b) 

r 

Given the energy levels and the total particle number the 
chemical potential is determined at any temperature and reso- 
nance energy by the constraint (O. 
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Figure 1 : (Color online) The density of states of a 3D optical lattice calculated by binning the energy levels into intervals of length A_B 
O.lSER.a- The blue (dark gray) and red (light gray) curves show the density of the atomic and molecular states, respectively. 



A. Energy levels in the lattice 

We consider a simple cubic optical lattice created by the 
overlap of three orthogonal pairs of counterpropagating lasers 
with wavelength Xl — 27r/fcL- This gives rise to the potential 



K(x) = l/o,.L(x), 



(9) 



for the particle type a (corresponding to atoms "a" or 
molecules "m"), where L is the dimensionless shape of the 
potential 



L(x) ~ [sm^{kLx) + sin^(fciy) + sm^{kLz)] 



(10) 



The molecules experience a lattice potential twice as deep as 
that in which the atoms move, Vb.m = 2Vb,a = 2Vb, since the 
Stark shift of a molecule is the sum of those for each atom. 

Now the problem is to understand how the difference in par- 
ticle mass and apparent lattice depth for atoms and molecules 
affects the behavior of their respective single particle states, 
obtained by solving the time-independent Schrodinger equa- 
tion 



2m„ 



+ Fo,.L(x) C(x)=i?;C(x). 



(11) 



There is a simple relationship between the atomic and molec- 
ular spectra, which is revealed by transforming ( fTTT ) for 
each species into their respective recoil energies, i?_R,cr = 



Vo,.L{^)\VA^) = Kfa{^). (12) 



where barred quantities are in recoil units. The advantage of 
Eq. ( fT2] i is that, since the left hand side operator only depends 
on a via the quantity Vb,o-, the spectrum is of the form 



e: 



(13) 



The relationship between Vb.a and Vb.m is quite simple: the 
molecules are twice as heavy as the atoms and see a lattice 
that is twice as deep as the one experienced by the atoms. As 
a result the atomic and molecular potentials are related as 



Vr 



0,m 



0,a, 



(14) 



i.e. in the respective recoil units the molecules see a lattice 
that is four times deeper. In general, an optical lattice with 
depth exceeding one recoil unit has a considerable effect on 
the spectral properties of the confined particles. Therefore, 
the difference in lattice depth for the atoms and the molecules 
has a rather profound effect on the properties of our system. 

We note that our choice to define the atomic energy origin 
as the lattice potential zero, rather than the atomic ground- 
state, is to emphasize the differential lattice confinement ef- 
fects on the atoms and the molecules. 



III. NUMERICS 

Because we are considering a separable potential ( fTOl i. the 
3D spectrum of the time-independent Schrodinger equation 
is most efficiently calculated via the ID eigenvalues, which 
are easily determined by numerical diagonalization. For the 
data presented here we have used a lattice with 31x31x31 
sites; this is sufficiently large that this system can be regarded 
as being approximately in the thermodynamic limit. We only 
exemplify results for f] = ^ since these represent the generic 
behaviour of the systems with a filling fraction less than unity. 

We numerically determine the chemical potential as a func- 
tion of T and E^cs under the condition of conservation of the 
total number of particles, (O. For a given point (/i, Eycs,T) 
in the phase diagram any thermodynamic quantity, such as the 
energy or the entropy, may then be calculated as the sum of an 
atomic and a molecular contribution. 
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A. Density of states 



IV. PHASE DIAGRAM 



The exact 3D energy levels can be found by making all pos- 
sible combinations of three ID energy levels. However, for the 
purposes of calculation the number of individual states needed 
is unwieldy, and it is desirable to construct a binned density of 
states by gathering the energy levels in small energy intervals 
and representing them by the centre of the interval. Doing 
this we obtain the graphs shown in Figure [T] The reliability of 
the binning procedure is confirmed in Figure \lla.), where we 
recover the usual dependence of the densities of states 
in the limit of vanishing lattice depth. By comparing Figures 
[nb) and (d) one sees that the molecular density of states for 
Vq = 2£'/f a and the atomic density of states for Vq = SEu^a 
have the exact same shape, though the latter is scaled on both 
axes by a factor of two (see Eqs. (flJl i and (fT4]i). To simpli- 
fy the numerical task without losing effects due to details in 
the energy spectrum we have used the exact energy levels for 
the lowest energy bands only and the binned density of states 
for the higher energy levels. The justification for this approx- 
imation is that the high energy domain is only relevant if the 
temperature is high, in which case the distributions ,fa{E) are 
sufficiently slowly varying that we can consider them constant 
over small energy intervals. 

As the lattice depth is increased, band gaps emerge in the 
3D density of states due to the gaps in the ID energy spec- 
trum. With reference to results shown in Figure [T] we make 
the following observations of the spectral properties: 



The molecular density of states has clear characteris- 
tics of a being in a much deeper lattice than the atom- 
ic system, i.e. for a given Vq the molecular spectrum 
has smaller band widths and larger band gaps than the 
atomic system. Most importantly, even for shallow lat- 
tices (e.g. Vq = 2Ejf^a), the molecule ground band 
is very narrow compared to the atomic ground band. 
These observations are consistent with the discussion 
below Eq. ( fT4] i. 



2. There is a positive offset in the ground state energy of 
the molecules relative to the atom states. As shown in 
Appendix lAl this shift arises from anharmonic effects 
in the lattice, and for the deep lattice limit it is given by 



The phase diagram is characterized by the molecule frac- 
tion, defined as 



E^„ 



8 xAW^ 



ERa- 



(15) 



Such shifts should be measurable in experiments as a 
displacement of the magnetic field position where half 
of the atoms have been converted to molecules in an adi- 
abatic sweep across the Feshbach resonance (see sec- 
tion[V]i. 



X{Eres,T,Vo,r]) 



2N„ 



(16) 



ranging from zero (all atoms) to unity (all molecules). In this 
section we first examine the molecule fraction as a function 
of the position of the resonance and the temperature of the 
system. We then consider how the low temperature chemical 
equilibrium is affected by changing the lattice depth. 

A convenient energy scale is the Fermi energy E^ir], Vb), 
which is taken to be the highest occupied energy level, when 
all the atoms are unbound. We note that for our choice of en- 
ergy origin the relevant Fermi temperature for characterizing 
degeneracy is given by Tp = {E-p — E^)/kB, with E^ the 
atomic ground state energy. 



A. Properties at fixed lattice deptli 

The molecular fraction is shown in Figure |2] as a function 
of E^cs and T for a fixed lattice depth. Those results show 
that tuning the Feshbach resonance, i.e. i?ios, provides a di- 
rect way of varying the composition of the gas. For large and 
positive values of E^-cs, the molecular states are at much high- 
er energy than the atomic states and are not thermally accessi- 
ble, thus realizing a pure atomic gas. Conversely, for large and 
negative values of Eyes, the system exists as a pure molecular 
gas. At low temperature the transition between these two lim- 
iting regimes occurs when Ey^s is close to Ep (see Figure |2]i. 
We remark that for a free system or a gas trapped in a har- 
monic potential the transition from atoms to molecules starts 
at E,,s ^ 2Ef UM- 
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Figure 2: (Color online) Molecule fraction in a {Ercs,T) phase dia- 
gram for Vo = ^En^a and 77 = i. 
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Figure 3: (Color online) Schematic illustration of how fi changes when the resonance energy is varied across Ep- For clarity we consider 
a deep lattice where the density of states of the lowest atomic- and molecular band is a narrow peak as illustrated by the solid lines. The 
distributions fa{E) and fm{E) are indicated by the dashed lines, and the shading under the peaks indicates the population of the atomic 
and molecular levels. Note that the iSrcs-dependence is included in the position of the peak in the molecular density of states and not in the 
Bose-Einstein distribution. 



The characteristic asymmetric fan-shape of the x-contours 
occurs because when the temperature increases, molecules 
will dissociate as two seperate atoms are entroptically favor- 
able over a single molecule. Thus if we follow vertical lines 
in the phase diagram in Figure|2]towards higher temperatures, 
the fraction of molecules must decrease eventually. 

The chemical potentials for degenerate Bose and Fermi gas- 
es exhibit markedly different behavior, the former being con- 
strained by the lowest available single particle level, while the 
latter depends explicitly on the number of particles in the sys- 
tem. Furthermore, in our system chemical equilibrium con- 
strains the molecular chemical potential to be twice that of the 
atoms. Thus it is clear that the change in atomic and molecu- 
lar populations induced by varying the applied magnetic field 
must be accompanied by a change in the behavior of the chem- 
ical potential. This can be described qualitatively by consid- 
ering how the atoms and molecules are distributed over their 
respective energy levels when Eras is swept across the Fes- 
hbach resonance as sketched in Figure [3] We identify three 
separate regimes: 

(a) Pure Fermi gas limit, E^es oo: If E^-cs is sufficiently 
high, the lowest molecular energy level E^^ + E^-cs is 
nearly unoccupied and all atoms are unbound. We then 
have 



fj. : 



Ev for T « 0. 



(17) 



This corresponds to Figure[3]!a) where the narrow ener- 
gy bands are sketched as peaks. 



For filling fraction ?/ < 1 the Fermi energy lies in the 
ground band. In the limit of a deep lattice the lowest 
band becomes sufficiently narrow that we can make the 
approximation E-p « E^ (see AppendixlAl). 

(b) Intermediate region, Ey^s ~ E-p: When E^ + E^^s ap- 
proaches i^L,n from above, molecules will start to form 
while the number of atoms decreases as shown in Fig- 
ure [Sjb). When i?res decreases further, i^l^ also has to 
decrease due to the condition (|6]l. Since = 2/i, 
the atomic chemical potential therefore also starts to de- 



(c) Pure Bose gas limit, iJres ^ —oo: Finally, when the 
resonance energy becomes low enough, /i has decreased 
so much that almost no atomic states are occupied, and 
the result is a nearly pure molecular gas corresponding 
to Figure[3tc). In this limit we have « E'^^ + Ey^s 
at low temperatures, i.e.. 



for T « 0, 



(18) 



From this description it is clear that the transition from atoms 
to molecules is more abrupt the narrower the lowest energy 
bands are. 

The low temperature behaviour of /i at half filling is shown 
in Figure |4] for a range of different potential depths and the 
effects discussed above are evident. Linear fits to /i in the 
molecular regime turns out to be in excellent agreement with 
the pure Bose gas limit estimate (fTsT l and the /^-plateaus in 
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Figure 4: Chemical potential as a function of _Eres for r] = \ and 
different potential depths in the limit T 0. This corresponds to 
following a horizontal line in the bottom of phase diagrams like in 
Figure|2] The linear behaviour in the left side is characteristic for the 
pure Bose gas regime, and the difference in the slopes arises because 
the Fermi energy varies with the potential depth. Inset: chemical 
potential curves for Vo = 5, 10, 15, and 20i5K,a (from bottom to 
top) for a finite temperature of ^.QEn^a/ks- 



the pure Fermi gas limit reveal the dependence of Ey on the 
lattice depth. Note that the arguments above do not tell us any- 
thing about the behaviour in the transition zone between the 
constant and the linear regime, but merely that at sufficiently 
high E^cs the chemical potential must be constant, while for 
E^cs tuned sufficiently below the resonance the chemical po- 
tential is linear. The fact that the transitions between these two 
regimes is so sharp as in Figure|4]is a consequence of the very 
low temperature. However, the general behaviour of plateaus 
on the atomic side of the resonance and a linear variation of 
/I on the molecular side remain also at higher temperatures in 
the limits E^^s ^ ±cx3 as illustrated in the inset of Figure H] 



V. DISSOCIATION ENERGY 

A noticeable feature in Figure |4] is that the transition point 
where the atoms start to form molecules varies with the poten- 
tial depth. In the absence of the optical lattice potential it lies 
at E^cs ~ 2i?F, as is well-known from BCS-BEC cross-over 
theories for a narrow Feshbach resonance ifTTIl . As the lattice 
depth is increased the resonance energy corresponding to the 
onset of the transition first decreases and then increases slight- 
ly towards a limiting value of E^^s ~ E-p in the deep lattice 
limit. We will elaborate on these effects in the following and 
explain the behavior for deep lattice potentials. 

To quantify the location of the chemical transition we de- 
fine the dissociation energy i^dis as the value of the reso- 
nance energy where half of the atoms have been converted 
to molecules, xC-E-dis) = 5 lUl- In addition, we introduce the 



zone, such that no unbound atoms are found for E'les < 
and no molecules exist for E^^s > E'^^^. The lower limit is 
defined by the condition 

e2„ 



e: 



2Et 



(19) 



which specifies where the first molecules start to break up as 
the resonance energy is increased from the molecular side of 
the resonance. The other end of the transition zone, where 
the first molecules are formed as the resonance energy is de- 
creased starting with a pure atomic gas, is defined by the con- 
dition that the lowest molecular level passes twice the energy 
of the highest occupied atomic state, i.e. 



(20) 



A. Free space limit 

These estimates simplify in the limit Vq = where the 
lowest atomic and molecular energy levels are only separated 
by E,^^, i.e. E^ = 0. We then have 





2Ey 



for 



= 0. 



(21) 



To find i?dis for a vanishing potential we note that in this case 
the density of states has a squareroot dependence on the ener- 
gy 



Pa{E) = CVE 



(22) 



with C = ttAP /2E'^^.The following argument applies at 
zero temperature: the number of unbound fermionic atoms, 
found by integrating the atomic density of states from zero to 
the chemical potential, is Na ~ "^Cp?^'^ . \i p — Ep, this 
equals the total number of atoms in the system, and in general 
we obtain the expression 



\ 3/2 

Ef) ' 



(23) 



for the molecule fraction (at T = 0). In the transition zone 
Ercs ~ 2[i (since E%^ ~ 0), such that 



X = 1 



/ Er, 



\2Ef 



3/2 



(24) 



and the dissociation energy in this limit then follows from the 
condition x ^ 



Edis = 2^/^Ef ~ 1.26£;f for Vn = 0. 



(25) 



This is confirmed by the numerical calculations, see Figure|5] 
Similar arguments can be used to show that ii^dis = 2^/'^£'f 
for a gas trapped in a harmonic potential j^). For a harmonic 



quantities iJ^g, delineating the zero temperature conversion 



trap one also finds E^-^_, 
El « i?F. 



and E 



dis 



2Ev since E, 



•0 




B. Deep lattice limit 



In the opposite limit of a deep optical lattice the 3D ener- 
gy bands become increasingly narrow and E^^^ and E'^^^ ap- 
proach each other. It is therefore reasonable to estimate E'dis 
to lie halfway between them, and inserting ( fT9] l and (|20] | yields 



Edis Ef - - E° ) 



for 



Vo 



oo. 



(26) 
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In the deep lattice limit the narrow ground band is well- 
described by a tight-binding approach (see Appendix [All, 
which gives an analytic expression for i?,'], — E^ (IA13b . We 
then have that 



0.5 



Edis ~ Ep 



1 



ERn. 



(27) 



R.a 



This estimate, valid for a deep lattice, is plotted in Figure |5] 
where we also show the numerical halfway mark for the atom- 
molecule conversion, i?dis, and there is good agreement be- 
tween the numerical results and ( |27] | in the deep lattice limit. 
Note that since E-p « E^ for a deep lattice, (IA121 i indicates 
that in this limit Ep / E^^a ~ 3 ^/VoJEr^ and hence that E^is 
eventually approaches Ep as the lattice depth increases 



Ed 



Ep for Vb — > oo, 



(28) 



since the first term in dZTb is then the dominant one. The de- 
viation from this limit in a deep lattice is given by the anhar- 
monic corrections in ( |27] |. 



C. Lattice induced resonance shift 



Figure 5: (Color online) Dissociation energy vs. optical lattice depth. 
For deep lattice potentials Edis can be approximated by Eq. ilH 
(white line), while the limiting behavior for a vanishing lattice depth 
is given by Eq. l |25t (dashed white line). Note that since the Fer- 
mi energy increases with increasing lattice depth, the scaling on the 
vertical axis depends on Vq. 



For ^Li in a lattice with wavelength Al 



ing A/i = 2/iB we find that B 



- ^dis 



1032 nm and tak- 
0.11 G for a 



lattice depth of Vb = l5ER^a assuming a filling fraction of 
?7 = i . This is comparable to the magnetic field width of the 
resonance AB = 0.23 G at resonance at 543.26 G iH. A 
similar expression has been found experimentally for the lat- 
tice induced shift of the magnetic position where the Feshbach 
molecule enters the continuum LISJ . 



These results demonstrate that the lattice induced energy 
shift between atomic and molecular degrees of freedom alters 
the dissociation condition for the Feshbach resonance. In an 
adiabatic sweep of the magnetic field across the resonance the 
effect of the lattice is to modify the molecule formation curve 
(X vs. B). In particular, the magnetic field corresponding to 
conversion of half of the atoms to molecules, B'^"^, is shifted 
with respect to its value in the absence of the optical lattice, 
B^i^, by an amount - B^i^, which in the deep 

lattice limit may be estimated from 



E. 



Vo 
dis 



E° 



J^p 



l.26E^ 



A/i 



(29) 



Here we have used that i^ios varies linearly with B with a 
slope A/i given by the magnetic moment difference between 
the open and the closed channels. The Fermi energy with no 

2/3 

lattice is Ep = (^) Eji a and for deep lattices we can 
approximate E^" w S\/VqEr^. Inserting this we get 



\ fJ-B J V amu / \ nm / 



6jj \2/3 



1.43 X lO^^G. (30) 



VI. CONCLUSION 

We have presented results for the thermodynamics of an 
atom-molecule mixture in an optical lattice potential under the 
condition of chemical equilibrium adjusted by a Feshbach res- 
onance which controls the energy offset between the atomic 
and molecular levels. 

The phase diagram has been determined and we have ana- 
lyzed the behavior of the chemical potential and the fraction of 
molecules with emphasis on the low-temperature regime. In 
particular we have identified the dissociation energy defined 
as the energy offset where 50% of the atoms have been dimer- 
ized into molecules. The chemical conversion takes place in a 
transistion zone of the magnetic field controlled resonance en- 
ergy, which narrows and shifts as the lattice depth is increased 
and the Bloch bands approach discrete energy levels. Further- 
more, the dissociation energy shifts with the lattice depth, as 
the molecules and the free atoms experience different lattice 
potentials, and hence their energy levels are displaced with 
respect to one another as the lattice depth is changed. In the 
deep lattice limit an analytic, precise expression for the center 
of the transistion zone was obtained by treating anharmonic 
corrections to the energy levels in the lattice wells perturba- 
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lively. 

Our results thus show that in an optical lattice potential the 
position and width of the Feshbach resonance, as indicated by 
the conversion of atoms into dimers and vice versa, depends 
on the lattice depth. We have shown that the lattice induced 
shift should be measurable for an atom with a narrow reso- 
nance such as ^Li. 

By focussing on the ideal gas case we have been able to 
clarify the effect of the uniform optical lattice potential on the 
atom-molecule equilibrium. The inclusion of interactions and 
an inhomogenous trap potential will be areas of future devel- 
opment to completely describe this system and make quan- 
tatative comparison with experiments possible. 
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Appendix A: THE DEEP LATTICE LIMIT 

When the optical lattice is very deep, we can make a tight- 
binding approximation and regard the potential as a collection 
of independent wells. The lowest energy levels are then 
those for one well with a degeneracy equal to the number of 
sites. A Taylor expansion of the ID potential (Vb,(T sin^ (kLx)) 
at a; = to 6th order gives 



V^'^'ix) « Vo^. ( (kLxf - ^{kLxf + ^{kLxf ] , (Al) 



so each lattice well can be described as a harmonic oscillator 
with the frequency 



with two anharmonic terms 



Hi = ^(kLX) , 

He = -^(kLx) , 



(A2) 

(A3) 
(A4) 



that we will treat with perturbation theory. The unpertubed 
energy levels for the harmonic oscillator are 

= (n +^\na;^= {2n + 1) /^V^, (A5) 



and the first order corrections 



(A6) 



can be calculated by expressing the factor {k^xY in terms of 
the ladder operators for the harmonic oscillator 



[kLxT = 



1/4 



+ a) 



We thus obtain the anharmonic first order energy shifts 



AS. 



n(l) 



E 



A pn(l) _ 



2n^ + 2n+l 



(A7) 



(A8) 



1 E 



36 yOVA^ 



'\n^ + 6ri^ + 8n + 3) , (A9) 



but we also need to consider the second order pertubation con- 
tributions from Hi to the n'th energy level 



ra^n 



{m\Hi\n) p 

cn cm 

<--(T.2 <--(T,2 



I (to| (a^ + d)^ \n) 



(AlO) 



since this contribution is of the same order as the first order 
pertubation contribution from Hq. In the general case the sum 
only gets contributions from 77i = n ± 2 and m = n ± 4, 
and in the simplest case n = 0, where the contributions from 

m = n — 4 and m = n — 2 vanish, we get 



. ^0(2) _ 7 Er 



(All) 



resulting in the approximation 




4 48 ^Vo,a/En,, 



ER,a (A12) 



for the lowest ID energy level. 

This derivation is valid for both atoms and molecules, and 
therefore the splitting between i?" and E'^j in three dimen- 
sions is given by - E^ = 'iiS^^ - i.e. 



iEl-EO)^l(l + i 



81 8 ^Vn,a/En,a 



EB,a (A13) 



where we have used the relation (fT4] i and that En^m = 
i5fl_a/2. The last term in ( |A13t can be made as small as de- 
sired by choosing a deep enough lattice potential, but the first 
term always remains resulting in an unavoidable energy differ- 
ence between the lowest energy band of the two components. 
As discussed in section |IV] this offset has observable conse- 
quences for the onset of molecule production in sweeps of the 
magnetic field across the resonance. 
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